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Abstract 

Wc collect and reanalyze about 200 GRB data of prompt-emission with known redshift observed until 
the end of 2009, and select 101 GRBs which were well observed to have good spectral parameters to de- 
termine the spectral peak energy (E p ), 1-second peak luminosity (L p ) and isotropic energy (E- lso ). Using 
our newly-constructed database with 101 GRBs, we first revise the E p -L p and E p -E iso correlations. The 
correlation coefficients of the revised correlations are 0.889 for 99 degree of freedom for the E p ~L p cor- 
relation and 0.867 for 96 degree of freedom for the E p -E- lso correlation. These values correspond to the 
chance probability of 2.18 x 10~ 35 and 4.27 x 10~ 31 , respectively. It is a very important issue whether 
these tight correlations arc intrinsic property of GRBs or caused by some selection effect of observations. 
In this paper, we examine how the truncation of the detector sensitivity affects the correlations, and we 
conclude they are surely intrinsic properties of GRBs. Next we investigate origins of the dispersion of the 
correlations by studying their brightness and redshift dependence. Here the brightness (flux or fluence) 
dependence would be regarded as an estimator of the bias due to the detector threshold. We find a weak 
fluence-dependence in the E p -E- lsa correlations and a redshift dependence in the E p -L p correlation both 
with 2 a statistical level. These two effects may contribute to the dispersion of the correlations which is 
larger than the statistical uncertainty. Wc discuss a possible reason of these dependence and give a future 
prospect to improve the correlations. 
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1. Introduction 

There are several correlations between the rest-frame 
physical quantities of GRB and their luminosity (or 
isotropic equivalent energy Ei so ). The first report 
was a variability-luminosity correlation by Fenimore & 
Ramirez-Ruiz (2000) which states that more variable 
events are more luminous. The next one was the lag- 
luminosity correlation reported by Norris et al. (2000) 
and Schaefer et al. (2001). Each pulse in the prompt 
emission has a spectral time lag, that is, a time delay 
of the soft-band emission compared with the hard-band 
one. According to this correlation, the events with large 
spectral time lag are dimmer than ones with short time 
lag. These two correlations are based on the temporal 
behaviors of prompt emission. 

Several correlations concerning the spectral property 
have also been suggested. Lloyd et al. (2000) found a 
correlation between the spectral peak energy (E p ) and 
observed energy fluence. They also mentioned the pos- 
sibility of correlation in the rest frame of GRBs. Amati 
et al. (2002) exactly mentioned a very tight correlation 
between the E p and the isotropic equivalent energy (Ei so ) 
in the GRB frame (see also Amati et al. 2006). This E p - 
E- lso correlation was confirmed and extended toward X-ray 



flashes by Sakamoto et al. (2004) and Lamb et al. (2004). 
Independently, Yonetoku et al. (2004) reported similar but 
tighter correlation between E p and the 1-second peak lu- 
minosity (L p ) called E p -L p correlation. Moreover, using 
the GRBs with measured opening half-angle, Ghirlanda 
et al. (2004) found that E p strongly correlates with the 
collimation-corrected gamma-ray energy (E 7 ). Firmani et 
al. (2006) reported the correlation among £" p -i?i so -To.45, 
where Xb.45 is the time spanned by the brightest 45 % of 
the total counts above the background. 

These correlations can be used as cosmological tools 
to investigate the physical environment of the early uni- 
verse. Yonetoku et al. (2004) used the E p -L p correla- 
tion as a redshift indicator for 689 GRB samples ob- 
served by BATSE without known redshift, and derived the 
GRB formation rate. Based on this GRB formation rate, 
Murakami et al. (2005) investigated the cosmic reioniza- 
tion epoch and the metal enrichment by the population-Ill 
stars. On the other hand, they are also useful tools for ex- 
tending the Hubble diagram to probe the cosmological ex- 
pansion history (Takahashi ct al. 2003; Oguri & Takahashi 
2006; Ghirlanda et al. 2006; Schaefer 2007). Kodama et al. 
(2008) calibrated the E p -L p correlation of nearby GRBs 
with the luminosity distance measured by the Type la su- 
pernovae (see also Liang et al. 2008; Cardone et al. 2009). 
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They succeeded in extending the cosmic distance ladder 
toward the redshift of z ~ 6, and estimated the energy 
density of dark matter and dark energy in high-redshift 
universe beyond z > 2. Furthermore, Tsutsui ct al. (2009) 
improved the Ep~E lso and the E p -L p correlation intro- 
ducing another parameter named luminosity time defined 
as Tl = Ei ao /Lp. Using this newly discovered E p ~L p ~Tl 
plane, they constrained the amount of dark matter and 
dark energy in the early universe more effectively. These 
papers well demonstrate the validity of these empirical 
correlations as cosmological tools. 

However, Nakar & Piran (2004) gave an argument 
against the presence of the E p -Ei SO correlation. They 
insisted that about 40% of 751 BATSE GRBs without 
known redshift do not satisfy the E p -E- lso correlation even 
if they assume any redshift, and also mentioned that clear 
outliers such as GRB 980425 and GRB 031203 exist. A 
similar argument against the E p -E- lso correlation has also 
been made by Band & Preece (2005). They concluded 
that 88 % of GRBs detected by BATSE cannot satisfy 
the E p -Ei so correlation because the observed -Ep-fluence 
ratios of these events exceed its maximum value around 
z ~ 4. If these arguments are true, we cannot use these 
correlations as cosmological tools. Butler et al. (2007), 
using the Bayesian approach to estimate E° hs for a lot of 
Swift events, indicated that dim events close to the detec- 
tor sensitivity would make large scatter on the Ep-E- lso 
and E p -L p correlations, and that there is a significant 
threshold effect (see also Lloyd et al. 2000). 

On the other hand, Ghirlanda et al. (2005a) gave a 
positive argument on the presence of the correlations. 
They used 442 bright BATSE GRBs with the pseudo red- 
shifts derived from the lag-luminosity correlation (Band, 
Norris & Bonncl 2004) and obtained the E p -E[ so corre- 
lation with the slightly different power-law index and the 
larger scatter than the original one. They found that the 
chance probability of the revised correlation is 2.1 x 10~ 65 . 
Similar conclusion has also been derived by Bosnjak et al. 
(2005). Ghirlanda et al. (2005b) also checked the valid- 
ity of the Ep-Lp correlation using 442 bright GRBs with 
the derived redshift and confirmed the correlation with 
the same power-law index within the error in the orig- 
inal Ep-Lp correlation by Yonetoku et al. (2004). They 
found that the chance probability of the E p -L p correlation 
is 1.6 x 10 -69 . The E p -L p and E p -E- mo correlations are 
independently tested by Suzaku-WAM and Fermi-GRB 
(Krimm et al. 2009; Amati et al. 2009). 

These correlations obviously have large dispersion com- 
pared with the statistical fluctuations. The origin of this 
data scatter is unknown, and should be revealed because 
the reliability and the accuracy of the GRB cosmology 
highly depend on the dispersion of correlations. This may 
be an intrinsic property of GRBs, or due to some instru- 
mental threshold effect. Recently Li (2007) and Basilakos 
& Perivolaropoulos (2008) tested the redshift evolution of 
the Ep-Ei so , Ep-Lp, E p -E~ 1 and other correlations and 
found no significant evolution, while the statistical errors 
are relatively large. Ghirlanda ct al. (2008) studied pos- 
sible instrumental selection effects on the E v and Fluence 



plane in the observer frame. In particular, they concen- 
trated on the trigger threshold (the minimum peak flux 
necessary to trigger a given GRB detector) and the spec- 
tral analysis threshold (the minimum fluence necessary to 
determine the value of E p ) . They showed these instrumen- 
tal selection effects do not dominate for bursts detected 
before the launch of the Swift satellite, while the spectral 
analysis threshold may be the dominant truncation effect 
of the Swift GRB sample. Nava et al. (2008) found that 
the E p of the fainter BATSE bursts is correlated with the 
fluence and flux with a correlation slope flatter than the 
one defined by the known redshift samples. They showed 
selection effects are not responsible for the correlations. 
About 6 % of these bursts are surely outliers of the E p - 
E- lso correlation, whereas there is only one sure outlier on 
Ep-Lp correlation. 

In this paper, using our newly-constructed database, 
we investigate whether the E p -L p and the Ep-E lso cor- 
relations represent the intrinsic property of GRBs or ar- 
tifacts of the truncation effect due to the detector sensi- 
tivity. Then, we study the origin of the dispersion of the 
correlations. In particular, we consider two possible ori- 
gins, threshold effect and redshift evolution. The former 
is that samples which arc detected marginally above the 
detector threshold may cause some systematic error. This 
should be tested by dividing the samples according to how 
larger than the threshold the flux/fluence is. However, the 
threshold here should not be the naive one because we are 
focusing on samples with the spectral parameters. Thus 
we have to use so called spectral threshold (Nava et al. 
2008) which needs detailed simulations. Instead, we sim- 
ply study flux/fluence dependence of the correlations as 
an estimator for the systematic error due to the threshold 
effect. This would be justified because the spectral thresh- 
olds for many detectors are roughly the same as shown in 
Nava et al. (2008). 

The structure of this paper is as follows. First we 
present a database of 101 GRBs with known redshift and 
wcll-dctcrmincd spectral parameters observed by the sev- 
eral independent missions in section 2. This database will 
be useful to estimate the intrinsic property of GRBs and 
we update the E p -L p and Ep-E lso correlations in section 
3. In section 4, we examine the possible systematic effects 
of these correlations such as the dependence of the corre- 
lations on the flux and fluence, and the redshift evolution. 
Finally wc will give discussion and summary in section 5. 

2. Data Selections and Analyses 

We collected and reanalyscd about 200 GRBs obtained 
by several instruments aboard the independent missions. 
Each instrument covers a different energy range with a 
different time resolution, so we should treat them carefully 
when we compare physical quantities such as E p and L p . 
In this section, we gather the observational data obtained 
until the end of 2009 and derive physical quantities by 
uniform criteria to construct a reliable database. After 
that, using this database, wc discuss the intrinsic property 
of the prompt emission of GRBs. 
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The prompt gamma-ray spectrum can be usually de- 
scribed as the spectral model of the exponentially- 
connected broken power-law function suggested by Band 
ct al. (1993): 



N(E) 



A 



100 kev) ex P( 

A ( E W(a 
^\ 100 kcV J \ 1C 



E 
E a 

0)E O 



for 



100 kcV 



o-/3 



cxp(/3 — a) for 



Here N(E) is in units of photons cm _2 s _1 kcV~ 1 . This 
function has four parameters, the low-energy photon index 
a, the high-energy photon index /3, the spectral break 
energy Eq and the normalization A. The peak energy can 
be derived as E° hs = (2 + a)Eo, which corresponds to the 
energy at the maximum flux in the vF v spectra. In this 
paper, we simply denote that E p — _E° bs (l + z) in the rest 
frame of GRB. 

According to BATSE observations, the average proper- 
ties of low- and high-energy spectral indices are a <~ — 1 
and (3 ~ —2.25, respectively (Preece et al. 1998; Preece 
et al. 2000). The recent Fermi-LAT observations sup- 
port that the high energy power-law index is consistent 
with the typical value of (3 = —2.25 beyond 1 GeV energy 
range (e.g. GRB 080916C, 081024B, 090323, 090428). 
Therefore, in this paper, if the high energy index has not 
been measured, we assume f3 = —2.25 as a fixed value when 
we estimate the bolometric flux and fluence. 

Here it should be noted that the spectra of some GRBs 
arc well fitted by cutoff power-law model. If some events 
show in fact the cutoff spectra, we will overestimate their 
bolometric flux and fluence when we fit the spectrum 
with the Band function (Shahmoradi & Nemiroff 2009). 
However, as far as we know, there is almost no posi- 
tive evidence about the existence of a clear cutoff in the 
prompt gamma-ray spectrum. Although many GRB spec- 
tra have been fitted with the cutoff power-law model, most 
of them are also well fitted by the Band function (see 
Kaneko et al. 2006 and Pal'shin et al. 2008 (GRB 080913), 
Sakamoto et al. 2009 (GRB 090516) and Pal'shin et al. 
2009 (GRB 090812) for more recent events). Therefore, 
at present, the fitting results by cutoff power-law model 
may tend to underestimate their total fluence and flux. 
Furthermore, as shown in Kaneko et al. (2006), it often 
happens that (3 cannot be determined by the data and 
the cutoff power-law model is sufficient to fit, if the peak 
energy is close to the high-energy end of the detector sen- 
sitivity, or if the event is so dim that the number of high 
energy photons is very small. Interestingly, simulations 
in Kaneko et al. (2006) showed that, if the signal-to-noise 
ratio is relatively low, a spectrum with the shape of the 
Band function is well fitted by cutoff power-law model. 
Thus, the assumption of the Band function for all events 
including dimmer GRBs would be reasonable at present. 
Our database is uniformly constructed in this sense. 

For the purpose of the complete coverage of database, 
we adopt three event selection criteria. (1) We selected 
GRBs with known redshift observed until the end of 2009. 
We have about 200 samples under this criterion from GCN 
Circular Archive (Barthelmy 1997) and GRBlog (Quimby 



et al. 2003). (2) We selected samples whose spectral pa- 
rameters are well measured. (3) We used samples whose 
total fluence and 1-second peak flux are reported. Under 
these three criteria we have 101 GRBs. For these samples, 
we can estimate the bolometric energy and the peak lu- 

^niifeity/^yS'axtending the spectral parameters toward the 
appropriate ehWgy band for each GRB with known red- 

^sM Q Tirm§§ 'paper, referring the previous works, we use 
1-10,000 kcV energy band in the rest frame of GRBs when 
we calculate the bolometric energy and peak luminosity. 
Then, we have to convert the observed fluence (S^bs) and 
the 1-second peak photon- flux (P Pj0 bs) within the energy 
range between E min and E max of each instrument, into 
the bolometric fluence (Sboi) an d the bolometric 1-second 
peak energy-flux (-F Pj bol) a s 
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Here, the integration is performed between the energy 
range of 1/(1 + z) kcV and 10,000/(1 + z) kcV for the Band 
function N(E). This is equivalent to the fc-correction. 
Then the bolometric isotropic energy (Ei SO ) and the 1- 
second peak luminosity (L p ) of 1-10,000 keV in the rest 
frame of GRB can be simply calculated as 
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(4) 
(5) 



Here, d^ is the luminosity distance calculated with the 
cosmological parameters of (O m ,fiA) = (0.3,0.7) and the 
Hubble parameter of H = 70 km s _1 Mpc _1 . 

It should be noted that we have very small number of 
GRBs satisfied with all three selection criteria if we use 
the information observed by a single instrument for each 
GRB. For example, the Swift/BAT is a powerful instru- 
ment to measure the low energy photon index a while 
it is difficult to determine the high energy index /3, and 
the E° hs value is not available for most GRBs. However 
the hard X-ray instruments, e.g. Konus, HXD-WAM and 
RHESSI strongly support Swift/BAT by measuring the 
higher-energy part of prompt spectra. Thus we frequently 
merge the spectral information reported by independent 
missions to obtain the entire shape of the spectrum. 

However, there are some GRBs observed by the inde- 
pendent satellites whose spectral parameters are incon- 
sistent with each other. In these cases, we basically ex- 
clude these samples from the table because we can not 
find any reasonable reason to choose one. However, there 
are two exceptions, GRB090423 and GRB090424. Fcrmi- 
GBM has a great advantage for the spectral measurement, 
especially for the E° hs determination, thanks to its wide 
band energy coverage. We chose the Fermi-GBM data 
rather than the Swift-BAT data for the two events, which 
have inconsistent E° hs values by the two detectors. Here 
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lve peaK pnoton nuxcs ror amcrenr tunc in- 
terval normalized by the one at 1024 msec. We used the 
lightcurve data with 64 msec time resolution observed by 
Konus. Peak photon flux with shorter time interval is sys- 
tematically brighter than the 1024 msec one. In this paper, 
we correct the time scale of peak fluxes if the reported value 
is in different time scale such as 64 msec or 256 msec. 

we would like to emphasize that E° ha values given by in- 
dependent detectors arc consistent with each other for a 
large majority of events within 1 a error. In those cases, 
we chose the one with a smaller error for them. 

Almost all satellite teams report the fluence (/S bs) and 
1-second peak photon flux (P Pi0 bs) in the energy range of 
their own instruments. However, some Konus reports did 
not include the 1-second peak flux but only 64 msec or 
256 msec one. Then using the lightcurve data published 
by the Konus team, we reanalyzed them and estimated 
the 1 second peak flux. The typical time resolution of the 
Konus/ Wind lightcurve is 64 msec. We performed the re- 
binning of 16 channel time bins into 1024 msec one with 
the running average method. Although we have 16 de- 
grees of freedom (d.o.f) to choose the start point when we 
calculate the 1 second peak flux, here, the 1 second peak 
flux is determined as the brightest one. This is a quan- 
titatively important modification which has been missed 
in the previous works. In figure 1, we show relative peak 
photon fluxes of Konus events for different time intervals 
(64, 128 ,256, 384, 512, 640, 768, 896 and 1024 msec) nor- 
malized by 1024 msec peak flux for each event. It can be 
seen that the 64 msec peak flux is systematically brighter 
than the 1 second one by 60-70 % level on average. 

In figure 2, we show the E p -L p distributions of the 
Konus data. The blue and the red filled squares are the 
peak luminosity measured with 64 msec and 1024 msec 
time scale in the observer frame, respectively. The black 
open squares are the data of 1 second peak luminosity 
observed by the other missions. The black, red and blue 
solid lines are the best fit power-law functions for the same 
colour data with the fixed index, respectively. We can find 
that the blue solid line is apart from the other two lines, 
and the normalization of the blue line is 1.70 times higher 
than the black one. On the other hand, the red line is al- 
most same as the black one with the factor of 0.96. This 
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Fig. 2. The E p -L p distributions of Konus data. The blue 
and the red solid squares are the peak luminosity measured 
with 64 msec and 1024 msec time scale in the observer frame, 
respectively. The black open squares are the data 1-second 
(or 1024 msec) peak luminosity observed by the other mis- 
sions. The black, red and blue solid lines are the best fit 
power-law functions for the same colour data with the fixed 
index, respectively. The blue data systematically distribute 
higher than the black line with the factor of 1.70. On the 
other hand, the red points are consistent with the black line. 

fact indicates that the 64 msec peak luminosity makes sys- 
tematic dispersion in the E p -L p correlation. Several past 
results, i.e. Nava et al. (2008); Ghirlanda et al. (2009) as 
well as our results by Kodama et al. (2008); Tsutsui et 
al. (2009), were argued about the selection effect and the 
GRB cosmology using the 64 msec peak luminosity for 
the Konus data. At present, since the number fraction of 
the Konus data is about 25 % of the entire 101 samples, 
the definition of the time scale should be treated correctly. 
Therefore the newly constructed database in this paper is 
more uniform compared with the past database, and ap- 
propriate to examine the origin of the data dispersion on 
the E p -L p correlation. 

The data are summarized in table 4 in appendix ??. 
Additionally, we also show the data excluded from the 
database in table 5 and 6 because of several reasons. For 
example, we know that so-called low luminosity GRBs 
and outliers are really exist. In this paper, we call the 
data "outlier" which locates outside of the 3 a confidence 
region of the E P ~L P and the E p -Ei SO correlations (see sec- 
tion 3). These groups are summarized in table 5. In ta- 
ble 6, we show several events with large ambiguity on 
the rcdshift and the spectral parameters. Several famous 
short GRBs are also listed in the same table as a reference, 
since their peak luminosity is determined by the millisec- 
ond time scale, and we can not convert them into the 
1-second peak luminosity. In the following sections 3, we 
discuss the E p -L p and E p -E- lso correlation. After that, in 
section 4, we examine a flux dependence and an redshift 
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evolution effect for these two correlations with table 4. 



3. Correlations 

In figure 3, we show the E p -L p (upper) and the E p - 
E lso (lower) correlations for all events listed in table 4 and 
table 5. We used the equations 4 and 5 when we estimate 
the isotropic energy and the 1-second peak luminosity, 
respectively. The solid black straight lines are the best- 
fit power-law functions for good data listed in table 4. 
The black curves around the straight lines represent 3 a 
statistical errors defined as, 
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where a a and <jb are errors on the normalization and 
the power-law index when we express the correlation as 
L p (or E iso ) = AEp 3 , respectively. The dotted lines are 3 o 
confidence regions including the data dispersion (system- 
atic error). 

In order to reduce the effect of samples which have un- 
expectedly large systematic error or which constitute dif- 
ferent (unknown) family of GRBs, we have to identify 
and remove outliers. We identified the outliers as follows. 
First, we use all samples except the two low- luminosity 
GRBs and obtain a tentative best-fit relation. Then sam- 
ples which are more than 3-sigma away from the tenta- 
tive relation are removed and a new (and still tentative) 
best-fit relation is obtained. Performing these procedures 
itcratively, we finally obtain the true best-fit relation and 
outliers. Thus, the selection of the outliers is not arbi- 
trary. As a result, we found 6 outliers, represented by the 
black points in figure 3. 

Our method is similar to those adopted in the analyses 
of Type la SNe (e.g. Kowalski et al. 2008) and Cepheid 
variables (e.g. Riess et al. 2009). It may be instructive to 
note that the number fractions of outliers are about the 
same for GRBs, Type la SNe and Cepheid. 

In figure 3, each color means the long GRBs (red 
points), high-redshift GRBs (z > 6, blue points) and short 
GRBs (green points), respectively. Recently, Levesque 
et al. (2010) suggested a definition of short-hard/long- 
soft GRBs with a statistical method. Their definition is 
roughly (1) events with TgQ S /(l + z) < 1.0 sec are short- 
hard GRBs. (2) events with T 9 ° b 7(l + z) > 2.0 sec are 
almost long-soft GRBs. (3) events with 1.0 < T 9 ° bs /(1 + 
z) < 2.0 sec can not be clearly defined but events with 
E° hs > 100 keV and shorter time scale (close to ~ 1 sec) 
are likely to be short-hard GRBs. Here, Tqq s is measured 
as the time duration during the 90 % of total observed 
photons have been detected while the Tqq s values highly 
depend on the energy band and the sensitivity of instru- 
ments. In this paper, we refer to their criteria while the 
original bimodal distribution of T§q measured by BATSE 
(Fishman et al. 1994) is seen in the observer frame. 

We estimated the best fit power-law function for the 
data listed in table 4: 
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Here we included not only the statistical errors but also 
the data dispersion (systematic error) with weighting fac- 
tor when we estimated the best fit values and errors 
in these correlations. The data dispersions are esti- 
mated as cr sy s,iogL p = 0.33 for the E p -L p correlation and 
0.37 for the E p -E- lso correlation, respectively 
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Then, the correlation coefficients are 0. 



for 99 d.o.f for 



the Ep-Lp correlation and 0.867 for 96 d.o.f for the E p - 
Ei SO correlation. These values correspond to the chance 
probability of 2.18 x 10~ 35 and 4.27 x 10~ 31 , respectively. 

The functional forms of these two correlations originally 
proposed in Yonetoku et al. (2004) and Amati et al. (2006) 
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respectively The power-law index of the revised E p -L p 
correlation is slightly different from the original one, but 
consistent within 2 a confidence level. Considering that 
Yonetoku et al. (2004) were able to use only 16 samples, 
2 a level agreement would not be strange. For the E p -E- lso 
correlation, the power-law index of our result is consistent 
with the original one by Amati et al. (2006) within 1 a 
confidence level. 

It is a hot topic whether the short GRBs are consistent 
with these correlations or not (Ghirlanda et al. 2009). In 
the criteria of short GRBs by Levesque et al. (2010), they 
are consistent with both correlations within the system- 
atic error. However the number of short GRBs is still 
poor, so more future observations and arguments should 
be required. 

4. Possible Origin of Data Dispersion 

As shown in the previous section, the two correlations 
are very tight, although there are some dispersion around 
the best-fit correlation: cr S ys,iogL p = 0.33 and cr S y S j og £ iBo = 
0.37, respectively. These correspond to the errors of about 
factor 2 when we estimate L p and i?; so from E p . 

One may wonder if these correlations represent the in- 
trinsic property of GRBs or just come from the trun- 
cation effect of the detector threshold (Butler et al. 
2007; Shahmoradi & Nemiroff 2009). In figure 4, we 
show the sample distributions in observer-frame quanti- 
ties (Ep hs -Fp^ ho i plane and E° hs -Shoi plane). We can see 
some correlations over three orders of magnitude in the 
observed brightness range, but these are substantially af- 
fected by the truncation of samples due to the detector 
sensitivity. In particular, GRBs with low brightness and 
large E° hs are generally hard to observe. Some authors 
insist that the truncation effect create the apparent cor- 
relations in GRB frame, such as the E p -L p and E p -E- lso 
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Fig. 3. The E P ~L P correlation (top) and the E p —E- iao corre- 
lation (bottom) for all events listed in table 4 and 5. The red 
and green points indicate the data of long-GRBs and short- 
-GRBs defined by Levesque et al. (2010) in the rest frame 
of GRBs. The 2 light-blue plots are the well-known low- 
luminosity GRBs of GRB 980425 and GRB 060218. The 
blue points are high-redshift GRBs with z > 6 (GRB 080913, 
GRB 090423). The solid black straight lines are the best fit 
function for long (red) and short (green) GRBs while two 
black curves around the straight line arc 3 a statistical er- 
ror. The dotted lines are 3 cr systematic errors in equa- 
tions 7 and 8. The 6 black plots are outliers which locates 
beyond 3 a confidence region from the best fit function of 
the E p -Lp and/or the E p -E iso correlations (GRB 050223, 
050803, 050904, 070714B, 090418, and 091003). 



correlations. 

However, this is not the case as we show. We divide the 
GRB events into three groups according to their bright- 
ness as shown in figure 4 and table 2 (for the details of 
the classification see the next subsection). As one can 
see, each group does not show significant correlation in 
the observer frame and it is evident that the truncation 
effect would be very small within each group. However, 
as figure 5 shows, each group shows a clear correlation in 
GRB frame. Thus, we can conclude that the correlations 
in GRB frame are not caused by the truncation effect. 

In this section, we seek for the origin of the dispersion 
of the correlations. The dispersion may be intrinsic and 
cannot be reduced, or some unknown systematic errors 
may contribute to the dispersion. Specifically we study 
the flux/fluence dependence and redshift evolution of the 
correlations and estimate the systematic errors accompa- 
nying them. 

4-1- Brightness Dependence 

Here we test the brightness dependence of the correla- 
tions using our database listed in table 4. The database 
consists of events with various brightness. The brightest 
event is GRB 991216 with the peak flux of F Pibo i = 5.80 x 
10 -5 erg cm _2 s _1 , while the dimmest one is GRB 020903 
with the peak flux of F Pi boi = 2.94 x 10~ 8 erg cm _2 s _1 . 
The difference is about 3 orders of magnitude. If the 
E p -L p and the E p -E lso correlations are affected by the 
detector sensitivity as suggested by Butler et al. (2007), 
systematic difference according to the brightness would be 
seen in the distribution of events on (E p ,L p ) and (E p ,Ei SO ) 
planes. 

We classify GRB events into three groups according 
to the bolometric peak flux to discuss the flux depen- 
dence of the Ep-Lp correlation, as shown in table 2: 
30 events with F p ^ Q \ < 5 x 10~ 7 (dim class), 41 events 
with 5 x 10~ 7 < -F p ,boi < 4 x 10~ 6 (middle class), and 
30 events with F p , bo i > 4 x 10~ 6 (bright class). We de- 
rive Ep-Lp correlation for each group and investigate 
whether they are consistent with each other. We per- 
form a similar analysis for the E p -E- mo correlation con- 
sidering three groups according to the bolometric fluence: 
27 events with Sboi < 5 x 10~ 6 (dim class), 41 events with 
5 x 10~ 5 < SW < 3 x 10~ 4 (middle class), 30 events with 
5*boi > 3 x 10~ 5 (bright class) as shown in table 2. 

In figure 5 (top left and right), we show the E p -L p 
and the Ep-E lso correlations for three brightness classes. 
The red, green and blue points represent the dim, middle 
and bright class as described in table 2, respectively. We 
adopt the power-law model for each group, and estimate 
the best fit function as well as the data dispersions of 
c S ys,iogL p an d Osys.iogBiso, respectively The solid line and 
curves for each colour mean the best-fit function and 3 a 
statistical boundary lines as same in figure 3. The results 
are summarized in table 1. 

In figure 5 (bottom left and right), we also show the 2 a 
acceptable regions for three groups. For the E p -L p corre- 
lation, we can recognize that three best fit lines are over- 
lapping. So we can say that the peak-flux dependence is 
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not significant in the E p -L p correlation. By contrast, for 
the Ep-Ei SO correlation, we found a weak trend that the 
events with larger fluence have the larger Ei so (or smaller 
E p ). 

The discrepancy between the bright and dim classes is 
over 2 a statistical level and would contribute to the dis- 
persion in the E p -Ei SO correlation. One should notice that 
the discrepancy is only a factor of 2 or 3, while the correla- 
tion itself covers beyond 3 or 4 orders of magnitude in the 
-Ei SO range. Therefore the detector sensitivity would con- 
tribute to the data dispersion of the E p -E- lso correlation, 
but not affect the argument on the existence of the corre- 
lation itself. We give an interpretation of this dependence 
in section 5. It should be noted that, as we pointed out in 
section 1, the flux/fluence dependence of the correlations 
can be regarded as an estimator for the systematic error 
due to the threshold effect. 

4-2. Redshift Dependence 

We also examine the redshift dependence of the Ep- 
Lp and the Ep-E lao correlations. This is a critical issue 
when we use these empirical correlations as cosmological 
tools like Type la supernovae. Our database covers a wide 
redshift range of 0.168 < z < 8.2 and we divide the sam- 
ples into three classes; 31 GRBs with z < 1, 36 GRBs 
with 1 < z < 2.5, and 35 GRBs with z > 2.5. In figure 6, 
we show the E p —L p (upper left) and the E p -E{ so (upper 
right) correlations of each class. The best-fit results and 
1 cr statistical errors are summarized in table 3. The 2 a 
confidence regions for three classes are also shown in the 
bottom left and right panels in figure 6, respectively. The 
meaning of each solid line is the same as one of figure 5. 

A difference between high- and low-redshift classes can 
be seen in the E p -L p correlation at 2 a level, and a sys- 
tematic redshift dependence is also seen in the power-law 
index. Especially the discrepancy is remarkable toward 
higher value of E p and L p . By contrast, in the E p -E- lso 
correlation, the three classes are consistent at 1 a level 
through the entire correlations. 



Fig. 4. The observed E p -Fj 



p,bol (left) and the Ep-Sbol 
(right) correlations. We classified the entire data into three 
groups according to each brightness range as listed in table 2. 
Each three group does not show the strong correlation be- 
tween E p and brightness, while the entire data set shows the 
clear correlation because of the truncation effect by the de- 
tector threshold. If each group consists of each E p -L p and 
Ep-Ei ao relation, then we may conclude the both correlations 
are real, and not due to the truncation effect. 



5. Discussion and Summary 

In this paper, combining data from multiple detectors 
and recalculating the peak flux for Konus events, we con- 
structed a uniform database of GRBs with known redshift 
and well-observed spectral quantities to determine E p , L p 
and Ei SQ . The 101 GRB samples enabled us to derive 
the Ep-Lp and E p -E- lso correlations with small statistical 
errors compared with their dispersions. By dividing the 
samples according to the observed flux/fluence and show- 
ing the correlation of each group, it was shown that the 
correlations are intrinsic to GRBs and not due to the trun- 
cation effect of the detector threshold. Then we examined 
the flux/fluence dependence and redshift evolution of the 
correlations, which is a crucial issue when we use GRBs as 
cosmological tools. We found a fluence-dependence in the 
Ep-E- mo correlation, and a redshift dependence in the Ep- 
Lp correlation. Both dependences are about 2 a statistical 
level and relatively weak so that they could still be used 
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Fig. 5. (Top) The left and right panels show the E p —L p cor- 
relation and the E p -E^ so correlation in the three flux (fluence) 
ranges listed in table 2. Red, green and blue points represent 
dim, middle and bright classes, respectively. The black line 
and curves are the best-fit function and 3 a statistical error 
derived in section 3 (see equations 7 and 8). (Bottom) The left 
and right panels show 2 a statistical error regions around the 
best fit functions of the E p —L p and E p —Ei ao correlations for 
the three flux (fluence) ranges listed in table 2. Each colour 
means the same of top panels. The black line and curves are 
the best-fit function and 2 a statistical error derived in sec- 
tion 3. The E p -L p correlation is consistent with each other 
in 1 a statistics while the E p —E- mo correlation show difference 
with 2 cr confidence level. 
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Fig. 6. (Top) The left and right panels show the E p -L p cor- 
relation and the E p —E^ so correlation in the three redshift 
ranges. Each colour represents z < 1 (red), 1 < z < 2.5 (green) 
and z > 2.5 (blue), respectively. The black line and curves 
are best-fit function and 3 rx statistical error derived in sec- 
tion 3. (Bottom) The left and right panels show 2 a statistical 
error regions around the best-fit functions of the E p -L p and 
E p -E iso correlations in the three redshift ranges. Each colour 
means the same of top panels. The black line and curves are 
best-fit function and 2 a statistical error derived in section 3. 
The E p —Ei 30 correlation is consistent with each other in f a 
statistics while the E p -L p correlation show difference with 
2 <r confidence level. 
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as cosmological tools. However, they would contribute to 
the data dispersion in the correlations. 

Let us give some comments on the work by Butler et 
al. (2007). They suggested, using the Swift data, the E p - 
E- lso correlation is due to the detector threshold. The E° hs 
value of many Swift events can not be determined by the 
observational data because of the narrow energy range of 
BAT instrument. So, using Bayesian statistics approach, 
they estimated the £"° bs values for 218 events based on 
a parent distribution of spectral parameters suggested by 
BATSE data (Preece et al. 2000). Butler et al. (2007) 
found a pscudo E p -E{ so correlation which is inconsistent 
with the original one by Amati et al. (2002). They claimed 
that this inconsistency is caused by the large intrinsic scat- 
ter of the correlation. In other words, the "tight" E p ~Ei so 
correlation is an artifact result, because almost all dimmer 
samples close to the detector sensitivity would be outliers 
of the E p -E- lso correlation by Amati et al. (2002). However 
we need caution to compare their work with others, be- 
cause their E° hs values are just estimated by the Bayesian 
statistics approach, and not observed. These simulated 
E"° bs values highly depend on the assumed parent distri- 
bution of spectral parameters. Our approach to prove the 
existence of the correlation is very different from this ap- 
proach. We do not assume anything about the statistical 
property of GRB samples. 

The fluence dependence in the E p -E- mo correlation may 
be consistent with the conclusion by Nava et al. (2008), 
although the statistical significance is only 2 a level. This 
dependence might be interpreted as follows. We can ob- 
serve only bright parts in the lightcurve for dim GRBs 
and tend to underestimate their total energies. In other 
words, the dispersion of the E p -E- lso correlation might be 
influenced by the instrumental threshold effect because 
the dimmer events can be observed only for the more 
sensitive instruments. In contrast, this effect would not 
be expected for the E p -L p correlation because it involves 
the peak luminosity (L p ) which is determined only by the 
brightest part of each event. 

We also found a weak redshift dependence in the Ep- 
Lp correlation. Of course, one interprets that this depen- 
dence is the intrinsic property of GRBs. However, let us 
point out a possibility of systematic overestimation of the 
peak luminosity L p for high-redshift GRBs. In the E p -L p 
correlation, we usually use 1 second peak flux in the ob- 
server frame when we estimate L p . However this means 
that the time scale of peak luminosity L p in the GRB 
frame depends on the redshift because of the cosmological 
time dilation. In other words, the time scale in the GRB 
frame becomes shorter for higher- redshift GRBs. In gen- 
eral, we can expect the peak luminosity increases as the 
time scale becomes shorter as shown in figure 1. For ex- 
ample, considering a GRB at z = 4, the observed 1 second 
peak flux corresponds to the 200 msec one at the GRB 
frame. According to figure 1, in this case, we systemati- 
cally overestimate the peak luminosity about 30-50 % on 
average. 

Thus, difference in the time interval in the GRB frame 
would induce a systematic error in L p . This effect would 
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make high-redshift GRBs look systematically brighter 
than low-redshift GRBs and may result in the apparent 
redshift dependence suggested in figure 6. Therefore, the 
redshift dependence of the E p -L p correlation might be un- 
derstood as a selection bias due to the inappropriate def- 
inition of L p , rather than an intrinsic property of GRBs. 
If this argument is true, we may enable to reduce the sys- 
tematic errors of E p -L p correlation when we use the same 
time scale in each GRB frame. Then we will make the Ep- 
Lp correlation more accurate luminosity indicator. This 
will be presented elsewhere in near future. 
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Table 1. The best fit results of the E p -L p correlation and 
the E p -E iao correlation in different flux (fluence) range. 



Class 




Ep-Lp 




a sysJogL p 




E -E- 

^lSO 






a sys,\og E iso 


dim 

middle' 

bright 


L v = 
Lp = 
Lp = 


1Q ^±um [E ■ 


1.82±U.lb 
1.48±0.16 
1.53±0.17 


0.23 
0.35 
0.35 


E- 

-^lSO 

E'iso 

E- 

-^lSO 


= 1() b2.82±U.U8 [£:? 
= 10 52.96±0.05 [£:? 
= 10 S3.30±0.11^ 


',355. 
',355. 
.,355. 


l.t)3±U.16 
1.24±0.13 
1.29±0.25 


0.24 
0.29 
0.43 



Table 2. The definition of three different flux and fluence 
ranges. 



Class 


F p ,boi range 


5boi range 




(erg cm~ 2 s^ 1 ) 


(erg cm~ 2 ) 


dim 


i^p.boi < 5 X 10-'' 


5 b oi < 5 x KT b 


middle 


5 x 10~ 7 < F p>bo] < 4 x 10~ 6 


5 x 10" 6 < 5 bo i < 3 x lO" 5 


bright 


F P .boi>4x 10" 6 


5 bo i > 3 x 10~ 5 



Table 3. The best fit results of the E p -L p and E p —E- lso cor- 
relations for three redshift ranges. 



Class 




Ep-Lp 






&sys 








&sys 


Z<1 

Kz<2.6 
z> 2.6 


L p 
Lp 
Lp 


= 10 ba.lH±0.07[^ 
= 10 52.47±0.06j Sj 
= 10 52.58±0.06[^ 


i,355. 
i,355. 
i,355. 


1.33±0.13 
1.54±0.15 
1.78±0.17 


0.32 
0.30 
0.23 




= io 52 - 97±0 - 08 [£; P ,355 : 
= io 53 - 05±0 - 07 [£; P ,355 : 


1.74±U.15 
1.44±0.21 
1.70±0.20 


0.40 
0.43 
0.31 
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Table 4. Spectral parameters of 101 GRBs with known redshift, which is applied for appropriate fc-correction. 
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l the observer frame, 
f Integrated between 1-10,000 kcV in the GRB frame. 

(t Short GRB with T g ° fi bs /{1 + z) < 1 sec in the rest frame of GRB and/or E p > 100 kcV (Lcvcsquc ct al. 2010). Greer 
i\ Merginal short GRB with Tg^ s /(1 + z) < 2 sec in the rest frame of GRB but E p < 100 koV. Green points in Fig.l. 
t High-redshift GRB with the rodshift of z > 6. Blue points in Fig.l. 
b Tg^ s highly depends on the energy band width of each instruments. 



Table 5. Spectral parameters of 2 low luminosity GRBs and 6 outliers excluded from the database of Table 4. 



t p,bol 



GRB 



rcdshift 



— FTTTB- 

-0.16 
5+0.42 
-0.42 
+ 0.37 
0.37 

I-H^O'. 1 ! 1 ! 

i.oo+^* x 

+0.1 
0.1 
+0.09 
0.09 
+0.01 
0.01 



, + O.U!) 

^ -0.09 

-2.25 + 7,x 
, + /" 
-Si 



(keV) 
Si 9 + "-^ 

69. 1+ 34 
+ 34 
34 
+ 335 
151 
+0.3 
0.3 
+ 780.0 
380.0 
+ 530.0 
-164.0 
+ 23.6 
23.6 



7S+1-01, 
7S -0.95' 
1S + 0.17) 
-0.17' 
59+0-^) x 10- 
+ 0.53 , 
-0.55' 
5 + 0.59, 
& -0.64' 

3 ± 0-.9 2 l>xl0- 

8i 2 : 6)xio- 

5 + 0.50, 
a -0.39' 



le 



' (5.97+J 
' (1.9li 9 
(6.99+0 
■ (1.77+0 
(3.72 + 
(5.07+1 
(2 53 + 
(5.05 + 



S bol 

+T 



(1.03+0-- 

(1 .07+0.: 
(1.61+0-; 
(S.48+ 2 -; 

(4.27+1-; 
(1.45 + 0': 

(i.i 8 +o; ; 

(2.78+0-- 



(org) 
4 + 0- 2 l, 



,/(.!+*) 
(sec) 



980425t 

050223* 

050803* 

050904* 

060218+ 

070714B* 

090418* 

091003* 



0.0085 

0.5915 
0.422 
6.295 

0.0331 
0.92 
1.608 

0.8969 



-0.99 



-0.86 
-1.3 
-1.13 



-2.06^ 



-2.25^ 



-2.50 



+ 0.10 
-0.10 
+ fix 
fix 



-2.25 
-2.25 + ^ 
-2.64+0- 



97 
436 
4.9 
1120.0 
610.0 
486.2 



9 5 )X10- 
') X 10" 



x 10- 



.19' 



.19 - 



54 



) X 10" 
) X 10" 
) X 10" 
) X 10- 
) X 10" 



;>x: 
i)x: 



10° 
10 4 



!)x 



f Low luminosity GR.B. 
X Outlier which locates beyond 3 cr from the best fit functii 



(i.05-- : 

(9.32+0-; 
(1.09+0- 
(1.60+0-; 

4+0.': 
-o.: 



(1.04^ 



39.7 
14.5 
77.4 

30.8 



21.5 
11.1 



. of the Ej. — L n 



the E n -E u 



No. ] Origins of Dispersions of ^-Brightness Correlation in GRBs 13 



Table 6. Spectral parameters of GRBs which 
excluded from the database of Table 4. 



have some ambiguities in the redshift or spectral parameters. These samples are 
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